arXiv: 1507.01409v2 [hep-ph] 22 Sep 2015 


Eur. Phys. J. C manuscript No. 

(will be inserted by the editor) 


The mass spectra and decay properties of dimesonic states, 
using the Hellmann potential 

Ajay Kumar Rai, D. P. Rathaud 

^Department of Applied Physics, Sardar Vallabhbhai National Institute of Technology, Surat, Gujarat, India -395007 


Received: date / Accepted: date 


Abstract Mass spectra of the dimesonic (meson - an¬ 
timeson) molecular states are computed using the Hell¬ 
mann potential in variational approach, which consists 
of relativistic correction to kinetic energy term as well 
as to the potential energy term. For the study of molec¬ 
ular bound state system, the Hellmann potential of the 
form T/(r) = — ^ -b — is being used. The one pion 
exchange potential (OPEP) is also incorporated in the 
mass calculation. The digamma decay width and decay 
width of the dimesonic system are evaluated using the 
wave function. The experimental states such as /o(980), 
6i(f235), hi(I380), ao(1450), /o(1500), /^(1525), /2(15 
65), /ii(1595), 02(1700), /o(1710), /2(1810) are com¬ 
pared with dimesonic states. Many of these states (mass¬ 
es and their decay properties) are close to our theoret¬ 
ical predictions. 

Keywords potential model ■ Exotic mesons • decays 
of other mesons 

PACS 12.39.Pn,14.40.Rt,13.25.Jx 

1 Introductiou 

The spectroscopy of hadrons imparts information of ba¬ 
sic force of nature. The theory of Quantum Chromo 
Dynamics (QCD) is a tool to understand the strong 
force in which hadron spectroscopy plays a key role in 
non-perturbative and perturbative regime. The new ex¬ 
perimental developments at Belle, BES, CLEO, CDF, 
LHC, BABAR brought out enormous data and came 
with large numbers of surprises [T1[21[31IH[S] . Number of 
new states are observed which do not fit in the con¬ 
ventional qq scheme, mainly in meson sector. All these 

'‘e-mail: raiajayk@gmail.com 


exotic states which do not fit in qqq and qq scheme re¬ 
quire extra theoretical attention [Hirflf^fl^lllO] . 

Recently, in partial wave analysis (PWA) oi J/tp —>• 
77777 ; BES HI collaboration observed most promising 
candidates for glueballs or multiquark structure below 
2.5 GeV [TTlfT^fT^ . In ■n~P ^ ryTr^ry channel GAMS 
group at GERN reported state near 1400 MeV [mu]- 
BNLE- 852 collaboration reported another state near 
1600 MeV in ■n~p reaction having promising hybrid 
structure [raiTClUTlfTB] . 

Meson spectroscopy, mainly, the light meson spectroscopy 
gives the information of the nonperturbative regime of 
QGD. The basic difficulty to study the light meson 
spectrum is that mostly resonances do not come out as 
narrow, isolated peaks [5]. Moreover, their large decay 
width make them difficult to identify experimentally as 
well as to differentiate them as threshold or as inter¬ 
ference effects, rather than pure resonances. Especially 
for scalar mesons, they are very difficult to identify due 
to large decay width, overlaps between resonances and 
background. Even, the scalar mesons have the same 
quantum numbers as of vacuum. Thus, the understand¬ 
ing the properties of scalar mesons may help us to un¬ 
derstand the mechanism of symmetry breaking 
Various mesons in low-lying sector like /o(980), ao(980), 
/o(1500) etc. observed experimentally, but their specific 
structure and properties are not understood theoreti¬ 
cally as well as experimentally. For example, the struc¬ 
ture of the scalar states /o(1500) and /o(1710) have de¬ 
bated since long time 

[21301131]. The authors in |191I201[^ indicated them as 

mixed or rather pure glue state respectively. Whereas, 
the Ref. P51[^l^[^l?7ll^[^ indicated /o(1710) as 
a Vector-Vector molecular candidate and /o(1500) as a 
glueball. The status of these scalar states are still a open 
question. The study of these exotic states may give the 
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answer to the basic questions of hadron masses, quark 
confinement, relevant degrees of freedom and interac¬ 
tion between constituents of the same [6117] . 

In the present study, we have investigated the molec¬ 
ular structures in the light flavour meson regime. Our 
dimesonic system consist of meson and antimeson {qq — 
qq). The masses of these states would be less than the 
sum of the masses of two constituent mesons mm- 
These loosely bound dimesonic system are similar to 
the deuteron like (Proton-neutron) bound state sys¬ 
tem. This approximation had been taken previously in 
Ref. [33] and introduced as ’deusons’. Furthermore, the 
molecular like picture had been studied by Silvestre et. 
al. [33] as diquonia. 

The multiquark states were predicted and studied pre¬ 
viously in MIT Bag model and in nonrelativistic poten¬ 
tial models jiUKin]. The molecular like structure have 
been studied and proposed in various models like poten¬ 
tial models p?^[H!niTlll[HKll!I^I!?71IMll391l4n] . chiral SU(3) 
quark model approach m. gauge invariant model m. 
Bethe-Salpeter equation approach [43], QCD sum rules 
[44] , an effective field theory approach [45] and nonper- 
turbative chiral approach [33] . We are using the poten¬ 
tial model to predict the masses, digamma decay and 
decay width of dimesonic systems in variational scheme. 
The potential model study is successfully explain the 
heavy flavour sector of mesons as well as multiquark de¬ 
scription of meson-antimeson or tetraquark states. The 
potential model are used to calculate the masses [511^1 
ITnil33ll34ll35l[36] and decay properties of various mul¬ 
tiquark states. Recently, the potential model has em¬ 
ployed successfully to study various spectroscopic study 
of mesons Whereas, various authors 

have successfully used non-relativistic potential model 
for the study of the molecular or tetra-quark structure 

[531135115311541155] . 

In our model for study of dimesonic states, we are us¬ 
ing the Hellmann potential with One Pion Exchange 
Potential (OPEP). The Hellmann potential is the su¬ 
perposition of coulomb and Yukawa potential. The su¬ 
perposed potential of coulomb and Yukawa was stud¬ 
ied by Hellmann long time ago [5i|571[55]. There are 
several authors who applied this potential for the cal¬ 
culation of the bound state systems for different values 
of their strength parameters [531I3U1I31 1162 1133] . All the 
phenomenological potential model approaches used by 
the authors, to study the mesons or multiquark spec¬ 
tra have dealt with complex interactions like, coulom- 
bic, confinement, instanton induced interaction, color, 
flavour and spin dependency to accomplished the bound 
state properties. At very short distances, the bound 
state system is produce due to delicate attractive and 


repulsive interaction. The Hellmann potential is sim¬ 
ply accomplished the complicated theoretical calcula¬ 
tion for the cancellation of the attractive and repulsive 
interaction at short distances, causes the small bind¬ 
ing energy of the bound state. The Hellmann potential 
with the OPEP used in the present model take care of 
short and long range behaviour of the interactions. As 
we aware that mesons are color neutral so we have not 
introduced the confinement potential for the dimesonic 
systems. It is well studied that the OPEP is mainly 
responsible for long range part of the strong nuclear 
force. The concept of meson exchange interaction has 
been used since long time for the N-N interaction [33] 
135] as well as for hadronic molecules [331I331I371I38] . 

To test the potential and mechanism employed, we first 
apply it to study the state /o(980) which was found 
in literature to be interpreted as molecular state m 
I38| . Moreover, the molecular like picture has been sup¬ 
ported by lattice QCD m- We calculate the spec¬ 
tra of /o(980) and obtained it’s mass, binding energy 
and root mean square radius. The results shows fairly 
good agreement with experimental results. Thus we are 
convinced and apply the same methodology to other 
dimesonic molecular states. 

The study of light meson sector needs relativistic treat¬ 
ment. We are incorporating relativistic correction to the 
kinetic energy term as well as to the meson-antimeson 
potential. The hydrogen like trial wave function is be¬ 
ing used for this study. We have calculated the digamma 
decay width and decay width of the dimesonic systems. 
The study of the digamma decay of scalars helps to 
distinguish among different scenarios for scalar meson 
structure m- The digamma decay of the dimesonic 
systems are calculated using the wave function at the 
origin [37]. For decay calculation, we adopt same for¬ 
mulation as used by Ref. [35], which was used to study 
the molecular picture of 1 ^ exotic states. 

The article is organized as follows. After the brief 
introduction, we present the theoretical framework for 
the study of dimesonic molecular states in semirelativis- 
tic approach with the Hellmann potential and OPE po¬ 
tential in section-!!. The calculation of digamma decay 
and decay width are discussed in section-HI. Then we 
presents the calculated results in section-IV and the last 
section is dedicated for conclusion. 

2 Theoretical Framework 

In the variational approach, we have solved the Schroed- 
inger equation by using trial wave function. The vari¬ 
ational parameter is being obtained by using the virial 
theorem. The (dimesonic) molecular system assumed to 
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Table 1 Experimentally observed states (in MeV) ^ 


States 

Mass 

Widths 

Possible 

Compared 


(MeV) 

(MeV) 


dimesonic 





states 

/o(980) 

990T20 

40-100 

0+(0++) 

K - K 

61(1235) 

1229T03 

142 ± 09 

1+(1+-) 

v-p 

/ii(1380) 

1386T19 

91±30 

o-(i+-) 

K - K* 

ao(1450) 

1474T19 

265T13 

l-(0++) 

p — u 

/o(1500) 

1505T06 

109T07 

0+(0++) 

p-p 

/'(1525) 

1525T05 

76±10 

0+(2++) 

p-p 

/2(1565) 

1562T13 

134T08 

0+(2++) 

OJ — U) 

/ii(1595) 

1594T15 

384T60 

o-(i+-) 

g' — u 

02(1700) 

1722 ± 16 

194T40 

l-(2++) 

K* - K* 

/o(1710) 

1722Aot 

135T07 

0+(0++) 

K* - K* 

/2(1810) 

1815T12 

197T22 

0+(2++) 

(f) — u> 


be consist meson-antimeson bound state. The Hamilto¬ 
nian of the dimesonic system is given by [7(T1[7T] 

H = ^ P'^ + ml + ^ P'^ + ml+ V{r) (1) 

where ma and mb are masses of mesons, P is the relative 
momentum of two mesons and V(r) is the molecular in¬ 
teraction potential of the dimesonic system. We expand 
the kinetic energy term of the Hamiltonian up to 0{P^) 


K.E. = 


pi 




/ 1 1 \ 
16 Vm® ml) 


/ 1 

+ 0{P^) 



( 2 ) 


V(r) is the molecular (meson-antimeson) potential, 

H (r) = (r)+(— + —'] {r)+o{ (3) 

\ma mbj \m^ J 

where. 


P(°)(r) = H(ri2) + H. 


(4) 


We consider the meson-antimeson interacting through 
the Hellmann potential. It is a superposition of the 
Coulomb and Yukawa potentials of the form [56ll57l[5^ 

O' P cri2 

V{ri2) = - -+B - (5) 

ri2 ri2 

Where Os and B are the strength of potential, C is 
the screening parameter and ri 2 is the relative coordi¬ 
nates between two meson and antimeson, the strength 
of coulomb potential is used as running coupling con¬ 
stant. While B is the Yukawa potential strength and it 
takes both positive and negative values. Ref. [60] had 
carried out a detailed study for different values of B and 
screening parameter C for low lying energy eigenvalues. 
For our dimesonic study, we have taken B positive. 


The value of the as running coupling constant is 
determined through the model, namely 


as{M^) = 


47r 


(11 - 


( 6 ) 


where M = 2ma mb/{ma+mb), Mb = 1 GeV, Aq = 
0.413 GeV and ny is number of flavour [72] . 

The nonperturbative form of VM ( 7 -) is not yet known, 
but leading order perturbation theory yields 


yd) (7.) = -CpCAallArl^- 


(7) 


where Cp =4/3 and Ca = 3 are the Casimir charges 
of the fundamental and adjoint representation, respec¬ 
tively m- The relativistic mass correction is found to 
be similar to the Coulombic term of the static potential 
when applied to the charmonium and to be one-forth 
of the Coulombic term for bottomonium |73j . Recently 
we have used this correction successfully to the study 
of Be meson m- 

The long range one pion exchange potential (OPEP) 
is used based on the assumption that molecular like 
structure of multiquark system is being deuterium like 
structure of nucleon P^I 66 ] . The OPEP for the NN- 
interaction takes the form [74] 


V 


Tr{NN) 


g^8 rn^ 
Att 3 


(n • o) 


Tn{r)Si2 + 





( 8 ) 


where gs is a pion-nucleon coupling constant, r and a 
are isospin, spin factors respectively. whereT^(r), V,r(r) 
and 512 are 


TAr) 



3 3 \ 

-'-TT ) 

m-^r m-^r‘^ J 




m-j^r 


(9) 


YAA 


r 


m-j^r 


( 10 ) 
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5'i 2 = 3(Ti ■ fcTj ■ rjr^ — <7i ■ aj 


( 11 ) 


The exchange meson having its own internal structure, 
to accomplished the finite size effect m, the usual form 
factor is appear due to the dressing of quarks and is 
assumed to be proportional to the exchange meson mass 

E3- 


TaAt) 


3 3 \ 

A^r ) A^r 


( 12 ) 


(r) 



(13) 


The functions T{r) and Y{r) with the finite size effects 
are takes the form 


Y{r) 

T(r) 


.1! 


Ai — m' 

7T TT L 


Y,{r) - Ay^Jr) 


mt 


via 


yl2 - m2 


T^[r)-^TAM 

m'j 


(14) 

(15) 


The more detail derivation of the potential can found 
in the Ref. m- The Si 2 is usual tensor operator, plays 
prominent role in NN-interaction. The matrix element 
of the tensor operator for L=0 vanishes [68]. If two 
hadrons are in an L=0 (ground state) state, the term 
with tensor operator vanishes as in our dimesonic case. 
Thus, the OPEP takes the form 


V 


TT 


l g^8 / rn^ \ 
3 47r \4:maTnt J 


{T^ ■ Tj) {at ■ aj) 



(16) 


where gs = 0.69 is a meson pion coupling constant, 
At^ = kuiTT, whereas, nia and rub are constituents masses, 
m-TT = 0.134 GeV and k = 2.2. For PV states the values 
of spin-isospin factor have been taken as (r^ • tj) (ui ■ Oj) 
= -3,1 for I = 0,1. While for VV states values taken as 
(■^i ■ ■^j) ('^i ■ ^j) = -6,-3,3 for isospin 1 = 0 and spin S 
= 0,1,2 respectively whereas for isospin 1=1 and spin 
S = 0,1,2 it takes values (r^ • Tj) {at ■ aj) = 2,1,-1. The 
values of spin-isospin factor are from Ref. [33] and the 
other parameters used in pion exchange potential are 

from Ref. [ 5 S 1 ISZ 1 ISH]- 

For the hyperfine splitting, the spin dependent in¬ 
teraction potential added perturbatively and takes the 
form 

VsD = l^^Si-S2\mf ( 17 ) 

9 marrib 

The spin factor 81-82 can be found by general formula 
Si-82 = i[(8i + 82)2-8i^-82^] m 


We have used the hydoginic like trial wave function 
such as. 


Rniir) 


f -1 - ly. y 

\ 2n{n + l)^\ ) 


e 2 




where fi is the variational parameter and 
is the Laguerre polynomial. The ground state mass of 
the low-lying dimesonic states is calculated by obtain¬ 
ing the expectation value of the Hamiltonian and the 
variational parameter(^) is obtained for each state by 
using the virial theorem [Tomii. 


Htjj = Eip and 


(K.E.) 


1 / rdV (r) \ 

2 \ dr / 


(19) 


We have included the kinetic energy term upto P® and 
potential with correction (Eq(3)) for the determination 
of g,. We fix color screening parameter C = 0.35 GeV 
for all combination and B = 0.4 for PP and PV-states 
and B = 1.2 for VV-states while strength of coulomb 
potential a(M^) is calculated as per Eq(6). The exper¬ 
imental (PDG) masses of the mesons are used for the 
present study, tabulated in Table-(II) [I]. 

The angular momentum,parity, spin, isospin, all are 
conserved in strong interaction, being good quantum 
number for dimesonic system. The orbital angular mo¬ 
mentum of meson and antimeson are Li and L 2 respec¬ 
tively. In same way spin denote by Si, S 2 and isospin 
Ii, l 2 - Employing the coupling rules, one have the rel¬ 
ative orbital momentum and total spin of system L12 
and S 12 whereas the total angular momentum to be J 
= L12+ S12 whereas I = Ii-I- I2. The charge conju¬ 
gation of the two particle system (meson-antimeson) is 
given by C = (_i)-£'i2-i-Si2 ijg the parity be P = 

Pi P 2 (—1)'^^^ [M] whereas G-parity is defined as G = 
(—The quantum numbers for conventional 
mesons as well as for exotic states have been discussed 
in Ref. [MUTS] . 

The one pion exchange potential (OPEP) for dimesonic 
system in Eq.(8) is spin-isospin dependent. Tdrnqvist 
in [33] had discussed the spin-isospin dependency of 
OPEP for deusons. The spin-isospin factor gives in¬ 
fluence to the potential and decides whether channel 
becomes attractive or repulsive. For VV-states we got 
repulsive channel for (S,I) = (0,1)(1,1)(2,0) while at¬ 
tractive channel for (S,I) = (0,0)(1,0)(2,1). Whereas 
in the case of PV-states the channel becomes attrac¬ 
tive for (S,I) = (1,0) while repulsive for (S,I) = (1,1). 
Due to parity violation, the two Pseudoscalar could not 
be bound by a Pseudoscalar 1331133] . Thus, we could 
not consider OPEP for PP-states. For such PP-states, 
the Hellmann potential given in Eq.(5) is being used 
with relativistic correction. The calculated masses are 
close to experimental measurements(PDG), tabulated 
in Table-(III). 
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Table 2 Masses of mesons (in MeV) ^ 


Meson rj rj' p u) K* ij) 

Mass 493.6 497.6 547.8 957.7 775.4 782.6 895.9 1019.4 


We have used Pseudoscalar and Vactor meson for 
the dimesonic system, having possible combinations (i) 
Pseudoscalar-Pseudoscalar(PP-state) (ii) Pseudoscalar- 
Vector (PV-state) (iii) Vector-Vector (VV-state). 


3 Di-gamma width and decay width 


TUq and q are the mass and 3-momentum of ex¬ 
change mesons (for present work, we consider only a 
pion). A is the free parameter which care of the off-shell 
effects at the vertices because of the internal structure 
of the mesons. The digamma decay widths are tabu¬ 
lated in Table-Ill along with masses. The decay width 
of compared states as well as our predicated states are 
tabulated in Table-IV. 


The digamma decay of the dimesonic molecules are esti¬ 
mated, using the wave function at the origin, in analogy 
to the two photon decay of parapositronium |37] . The 
di-gamma width for all mass combinations is given by 

(20) 

rriamb 

a = IAk indicate fine-structure constant, the factor ^ 
“"Tf’ masses of constituent mesons. 


The decay width is calculated using the wave func¬ 
tion at the origin. The dimesonic molecular state (con¬ 
stituent mesons a and b) is decays into mesons c and d. 
In the decay through meson exchange, the amplitude is 
proportional to the wave function square at the origin 
[3H]. The formula of |V'(0)|^ given by [55 1 I75 1 I77] 


|V’(0)|^ = 2M < V'I dP/dr I V’> 


( 21 ) 


with M as reduced mass of the system, '0(0) includes 
the effects of all ranges in r. 

The decay width is 


r = 


IV^(o)l^^ 

167rm3 


\M\^ 


( 22 ) 


Where rrim is the mass of dimesonic molecule and 1 
is the magnitude of 3-momentum of the decay product 
and is given by 


p ^ (mj + mt + mj) 

4m^ 

(m^mg -k -f rngm^) 

2TOm 


(23) 


TUa, mb are masses of constituent mesons and me, 
md are masses of product mesons, whereas A4 is the 
amplitude 


M = 


/ 


(24) 


4 Results and Discussion 

By using the Hellmann potential, OPE potential and 
relativistic correction, we have solved the Schroedinger 
equation to extract the masses, binding energy, decay 
widths and digamma decay widths for dimesonic states. 

In various combination of dimesonic states like PP, PV 
and VV- states, we have calculated the possible molec¬ 
ular like structure in light flavoured sector. Our calcu¬ 
lated results are tabulated in Table (III - IV). 

To describe molecular picture of dimesonic states, we 
have incorporated the short range and long range be¬ 
haviour of the potential. As the constituent of the molecule 
(meson) itself color neutral and with the small bind¬ 
ing energy and large radius (r in fm) in comparison to 
constituent size, we do not need to introduced the con¬ 
finement potential and so we have not dealt with funda¬ 
mental quark gluon condense state in the present study. 
Thus, we have not incorporated any mixing scheme of 
gluonia and quarkonia. Still, the results obtained in 
present study have reasonably good agreement with 
experimental measurements and explain the molecular 
picture well. 

We have expanded the kinetic energy term up to 0{P^). 

In the series expansion of the kinetic energy. For v << 
c, the effect of the higher order term of the momentum 
(n > 2) is negligible, even more, the higher order 
term has poor convergence. While the expansion term 
up to P* does not have a lower bound. So, the usable 
expansion to incorporate the relativistic effect is being 
up to P® [71]. We have dealt with the systems which 
have constituent masses below 2 GeV. For the lighter 
constituents, the effect of the momentum has small. As 
the V < c tends to u << c, the effect of the higher order 
terms contributes less than 1%. The effects of the higher 
order terms (up to P®) are very small, still, it has vari¬ 
ation with increasing masses and contribute {v —>■ c) 
to the net kinetic energy of the system. Therefore, it 


— mq 
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Table 3 Mass spectra, binding energy, root mean square radius and digamma decay width of dimesonic [qq — qq) systems 
with their values. 


System 


p 

R(0) 

V < > 

B.E. 

Mass 

Expt. ^ 


Exp. [I] 

State 



(GeV) 

(Gept) 

(fm) 

(MeV) 

(GeV) 

(GeV) 

KeV 

KeV 


PP-states 











K -K 

0+ (0++) 

0.219 

0.072 

11.63 

-20.87 

0.974 

0.990 

±0.020 

0.285 

0 

/o(980) 

v-v 

0+ (0++) 

0.265 

0.096 

09.63 

-26.24 

1.069 


0.414 



n-v' 

0+ (0++) 

0.409 

0.185 

06.24 

-46.75 

1.458 


1.526 



v' - v' 

0+ (0++) 

0.628 

0.352 

04.06 

-85.22 

1.830 


1.803 



PV-states 











v-p 

1+ (1+-) 

0.210 

0.068 

12.12 

-63.82 

1.259 

1.229 

±0.003 

0.208 


61(1235) 

rj — Lu 

0- (1+-) 

0.209 

0.067 

12.20 

-63.22 

1.267 


0.203 



K -K* 

0- (1+-) 

0.206 

0.066 

12.36 

-62.60 

1.330 

1.386 

±0.019 

0.237 


6,1 (1380) 

K-K* 

1+ (1+-) 

0.209 

0.067 

12.19 

-63.51 

1.330 


0.247 



rj- 4> 

0- (1+-) 

0.253 

0.090 

10.08 

-73.50 

1.493 


0.361 



v' -p 

1+ (1+-) 

0.351 

0.147 

07.26 

-94.25 

1.639 


0.316 



rj' — u! 

0- (1+-) 

0.350 

0.146 

07.30 

-93.45 

1.646 

1.594 

±0.015 

0.312 


61(1595) 

11 ' - 4‘ 

0- (1+-) 

0.429 

0.198 

05.95 

-107.4 

1.869 


0.574 




VV-states 


P 

- P 

0+ (0++) 

0.192 

0.059 

13.27 

-55.39 

1.489 

1.505 

±0.006 

0.079 


/o(1500) 

P 

-P 

0- (1+-) 

0.197 

0.062 

12.92 

-55.00 

1.492 


0.085 



P 

- P 

0+ (2++) 

0.194 

0.060 

13.15 

-54.10 

1.500 

1.525 

0.081 

0.081 

A(1525) 









±0.005 


±0.009 


P 

-P 

1- (0++) 

0.194 

0.060 

13.11 

-54.25 

1.490 


0.082 



P 

- P 

1+ (1+-) 

0.195 

0.061 

13.07 

-54.40 

1.493 


0.082 



P 

- P 

1- (2++) 

0.196 

0.061 

13.00 

-54.70 

1.499 


0.084 



UJ 

— ZJ 

0+ (0++) 

0.201 

0.064 

12.67 

-55.77 

1.502 


0.089 



U 

— ZJ 

0- (1+-) 

0.199 

0.063 

12.79 

-55.31 

1.506 


0.086 



UJ 

— ZJ 

0+ (2++) 

0.196 

0.061 

13.01 

-54.40 

1.514 

1.562 

0.082 

0.70 

/2(1565) 









±0.013 


±0.14 


K* 

-K* 

0+ (0++) 

0.231 

0.078 

11.05 

-59.34 

1.724 

1 '7O(y~\~0 .OOG 
-*-• '^^-0.005 

0.102 


/o(1710) 

K* 

- K* 

0- (1+-) 

0.229 

0.077 

11.14 

-58.89 

1.729 


0.100 



K* 

-K* 

0+ (2++) 

0.225 

0.075 

11.32 

-58.01 

1.730 


0.095 



K* 

- K* 

1- (0++) 

0.226 

0.076 

11.29 

-58.16 

1.726 


0.096 



K* 

-K* 

1+ (1+-) 

0.226 

0.076 

11.26 

-58.31 

1.730 


0.096 



K* 

-W 

1- (2++) 

0.228 

0.077 

11.20 

-58.60 

1.736 

1.722 

0.098 

0.30 

02(1700) 









±0.016 


±0.05 



-4> 

0+ (0++) 

0.255 

0.091 

10.00 

-60.91 

1.970 


0.106 



(p 

- 0 

0- (1+-) 

0.253 

0.090 

10.07 

-60.49 

1.974 


0.104 



(p 

- 0 

0+ (2++) 

0.250 

0.088 

10.22 

-59.68 

1.982 


0.100 



P 

— Zj 

1- (0++) 

0.196 

0.061 

13.04 

-54.40 

1.497 

1.474 

±0.019 

0.081 


ao(1450) 

P 

— Zj 

1+ (1+-) 

0.196 

0.061 

13.00 

-54.55 

1.500 


0.082 



P 

— Zj 

1- (2++) 

0.197 

0.062 

12.93 

-54.85 

1.506 


0.084 




(p - 

p 

1- (0++) 

0.222 

0.074 

11.47 

-57.82 

1.730 


0.070 

4, - 

p 

1+ (1+-) 

0.223 

0.074 

11.44 

-57.96 

1.733 


0.071 

4 - 

p 

1- (2++) 

0.224 

0.075 

11.38 

-58.25 

1.740 


0.072 

4- 

ZJ 

0+ (0++) 

0.228 

0.077 

11.17 

-59.10 

1.735 


0.076 

4 - 

ZJ 

0- (1+-) 

0.226 

0.076 

11.26 

-58.65 

1.739 


0.074 

4- 

ZJ 

0+ (2++) 

0.223 

0.074 

11.44 

-57.79 

1.747 

1.815 

±0.012 

0.071 


/ 2 ( 1810 ) 
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justified to incorporate the expansion up to P®. 

The correction used in the potential, have its dominant 
effect on potential energy. In our calculation, the effect 
of the correction to the potential energy part is increas¬ 
ing as mass of the system decreasing, and, it is about to 
be 25—40%. The effect of correction to potential has ap¬ 
proximate 38% in the lightest dimesonic system {K — K 
) of this work, while, the heaviest system (j)—4> has about 
28% contribution. We have calculated the partial decay 
width as per Eq.(22). The parameter A is assumed to 
be proportional to the mass of the dimesonic state, A = 
/C rrijn- We have chosen the Yukawa strength parame¬ 
ter and color screening parameter B = 0.4 and C = 0.1 
respectively to find mass spectra ior K — K molecule 
which believed to be /o(980). For the decay calculation 
oi K — K, we took 1C = 1.262. Our calculated results for 
mass, digamma and decay width are close to the PDG 
values [T]. Thus, we fix these parameters for all PP and 
PV-states calculation while for VV-state, we took B = 
1.2. For the dimesonic states having masses between 1 
GeV to 1.6 GeV, the constant /C = 1.085, whereas, for 
the states between 1.6 GeV to 2 GeV, 1C = 0.931. In 
PP-states, we used the Hellmann potential with rela¬ 
tivistic correction. Due to parity violation, the OPFP 
can not be apply to the PP-states calculation. 

^”i/o(980) = 0.974 GeV ; = 0.28 KeV ; 

-f^(/o(980)) = 94 MeV 

Our results for fo{980)fj^ is in good agreement with 
suggested mass and decay width (0.982±0.003 GeV and 
80±10 MeV) of Ref. [78]. The calculated digamma width 
of /o(980) is in close agreement with other theoretical 
predictions The TTTT decay mode of 

/o(980) is dominant and also agree with Ref [T]. For the 
molecular state, one should have small binding energy 
and large radius (compared to constituent energy and 
radius). Our results satisfied these conditions. The mass 
and decay properties are close to the experimental mea¬ 
surements. Thus, we are employing same methodology 
to calculate the properties of all the meson-antimeson 
combinations. The computed results are tabulated in 
Table-3 & 4. 

Pseudoscalar-Vector (PV) states :- 
In the PV-states, the ttV system being the lightest 
dimesonic state to be possible. In the present work, 
we have not considered ttV systems. The pion is too 
light as a constituent for dimesonic states which in¬ 
teract through exchange of pion itself. The pion as a 
constituent, carrying the large kinetic energy and it 
is difficult to overcome by potential energy to bound 
a molecule [33]. With the reasonably close mass spec¬ 


tra and quantum numbers, we have found three states 
5i(1235), /ii(1380) and /ii(1595) which could compared 
with rj—'p^K — K ,7]' —To dimesonic states respectively. 
The mass of 6i(1235)(^p) has fairly near to experimen¬ 
tal observation. 

^iTihi(i 235 ) = 1.259 GeV ; F^^ = 0.20 KeV ; 

^(&i(i235)) = 132 MeV 

the TTo; is experimentally observed dominant decay 
mode and our results are in good agreement. 

For /ii(1380) 

^mhi(i 380 ) = 1.330 GeV ; F^^ = 0.23 KeV ; 

l^(/ii(i380)) = 268 MeV 

For hi(1595) 

’7^m?,l(l595) = 1-646 GeV ; F^^ = 0.31 KeV ; 

-^(/ii(i 595 )) = 63 MeV 

In the case of the hi(1380)^^-j^j, our calculated de¬ 
cay widths are overestimated . Moreover, the mass of 
the state is underestimated around 50 MeV. For the 
/i(1595)(^/jj), the decay width is unacceptably under¬ 
estimated with experimental measurement with uit] ob¬ 
served decay mode. The mass of the state is also away 
from PDG values around 40 MeV. Thus, the obtained 
results indicates the rejection of molecular interpreta¬ 
tion of both the states. Thus, in the PV-states calcula¬ 
tions, the 6i(1235) is the only the state, has been found 
the strong candidature for dimesonic molecule. 

The Vector-Vector (VV) states :- 
In the VV combination, the dependency of one pion 
interaction potential on spin-isospin factor gives pos¬ 
sibilities of the numbers of combinations for dimesonic 
molecular like states mm- We have already discussed 
in the previous section about the matrix element, at¬ 
tractive and repulsive channels of the spin-isospin fac¬ 
tors. We have found bound states in repulsive spin- 
isospin channels, shows the dominance of hellmann po¬ 
tential and the mass correction over pion exchange. 
This can be clearly seen from binding energy tabu¬ 
lated in Table-IV for a state of same spin with dif¬ 
ferent isospin. With different spin-isospin combination, 
we have compared experimentally observed states with 
our dimesonic molecular states, with comparable mass 
spectra and quantum numbers. The states /o(1500), 
/^(1525), /2(1565), ao(1450), /o(1710), 02(1700), /2(18 
10) compared with different vector-vector combination, 
shown in Table-Ill. 

The /o(1500) is being compared as p — p molecule 
with (S,I) = (0,0), falling in 1500 GeV mass regime. The 
computed mass and decay widths are comparable with 
experimental values [33] and with Ref. [S3], suggested 
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mass 1.522 ± 0.005 GeV and decay width 108 ± 8 MeV. 
In Ref. [75] suggested decay width 131 ± 15 MeV, also 
in Ref. [8^l85ll86] . We have found the tttt decay mode 
is dominant. Experimentally two photon decay has not 
been seen yet. 

''T^m/o(i 5 oo) = 1.489 GeV ; = 0.079 KeV ; 

7^(?ii(i500)) = 110 MeV 

The state /2(1525) has been compared with p — p 
dimesonic state, with (S,I)=(2,0). Our calculated decay 
width for /2(1525) is consistent with Ref. [TIIST] whereas 
digamma width has agreement with Ref. [111571155] . 

"im/^(i 525 ) = 1-500 GeV ; = 0.081 KeV ; 

l^(/'(i 525 )) = 82 MeV 

The state /2(1565) is compared with oj—oJ molecule. 
The calculated mass of this state is underestimated 
around 40 MeV with the PDG value [I], and, the cal¬ 
culated digamma decay width is being underestimated 
about one ninth of the recent PDG value [T]. More even, 
the calculated decay width for observed decay modes 
are also far away from PDG values [T]. 

for /2(1565) 

Wm/ 2 (i 565 ) — 1.514 GeV ; F^-y = 0.082 KeV ; 

1"(/2(i 565)) = 81 MeV 

The states /o(1710) and 02(1700) are compared with 
K* — K dimesonic molecule with spin state S=0,2 re¬ 
spectively. The mass of /o(1710)^.-j^ has good agree¬ 
ment with PDG [T] and relatively near to Ref. [75] 
(1.750 ± 0.020 GeV) and Ref. [SnUSO] • Even, the com¬ 
puted decay width consistent with current experiment 
[T] and with Ref. [521111]. 

™iti/o(i 7 io) = 1.724 GeV ; F^^ = 0.10 KeV ; 

7"(/o(i 710)) = 143 MeV 

The mass of 02(1700)^.-^^ is in excellent agreement 
with experimental PDG value [T]. The digamma decay 
is found one third of the PDG value (0.30 ± O.QbKeV) 
[T], but, the decay width is in good agreement with ex¬ 
perimental values for observed decay modes [I] . 

777mo2(i7oo) = 1.736 GeV ; F^^ = 0.098 KeV ; 

-^( 02 ( 1700 )) = 181 MeV 

Further more, we compared p—ZJ and dimesonic 
molecules with the states ao(1450) and /2(1810) respec¬ 
tively. The mass of ao(1450)pjj is in agreement with 
PDG value [T]. Whereas the decay width of ao(1450)pjn 
has been found far away with PDG value [I], enforce us 


to rule out it as a dimesonic molecule. 

77imao(i450) = 1.49 7 GeV ; Fjj = 0.081 KeV ; 

7);ao(i450)) = 94 MeV 

Whereas, The mass of the /2(181O)0zj is off around 
50 MeV with experimental measurement [I], but the 
decay width are fairly comparable to PDG values in 
observed decay modes. 

^” 1 / 2 ( 1810 ) “ 1.747 GeV ; = 0.071 KeV ; 

-^(/ 2 (i 8 io)) = 196 MeV 

A.V. Anisovich et.al. [l^ made a remark on the sta¬ 
tus of the state / 2 ( 1810 ) such that the state may be 
confuse with /o(1790). [H] has made note, if there was 
not any confusion between /2(1810) and /o(1790) there 
should missing O^^ state. As we observed throughout 
our calculated results, all states (except few states) are 
underestimated approximate 20 to 40 MeV. If we con¬ 
sider it as limitation of our model, we can compare the 
state /o(I790) with (j) — oJ, with (S,I)=(0,0). So it may 
be possible that these two states /2(1810) and /o(1790) 
having molecular structure oi (jj — uJ with spin state S = 
0,2 respectively. But, the experimental status of both 
the states are not confirmed. Thus, it needs more atten¬ 
tion from theoretically as well experimentally. Whereas, 
our results favoured the molecular picture. 

The status of the scalar states /o(1500) and /o(1710) 
as a molecule, guleball or glue mixed state have de¬ 
bated since long time [Tmi^lTni^lMlIMlI^IMlITTlI^ 
lllllsisl]. The glueball and quarkonia with same quan¬ 
tum number and nearly close mass, may have interfered 
with each other and formed new state. If the one bare 
glueball and one or more quarkonia interfere, needs the 
mixing scheme to explain them. This interference or 
mixing of the state could be solved in mass matrix for¬ 
malism or in a linear combination of pure quarkonia and 
gluonia states by hnding the linear coefficients which 
mainly depends on decay properties |20] . The authors 
of Ref. [H] have predicated /o(1500) as a glue mix¬ 
ture with ss component and /o(1710) as predominant 
a pure glueball, also in Ref [HJUT]. In Ref. [1^ has indi¬ 
cated the /o(1710) as glueball candidate. Whereas the 
Ref. p51[Ml[^[^[?7l[^[51] have indicated /o(I7I0) as 
a Vector-Vector molecular candidate and /o(1500) as a 
glueball. So far, as per the present literatures, it has 
remained unclear and debatable the structure of the 
/o(1500) and /o(1710). Moreover, in the recent review 
of particle data [1] have noted on the status of the scalar 
states, in which /o(1710) together with /2(1525) were 
interpreted as bound systems of two vector mesons [T] 
m]. The molecular picture could be tested in radiative 
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J/'i/i decays as well as radiative decays of the states 
themselves OEOIISI]- While the vector component of 
/o(1710) might have the origin of the enhancement seen 
in near threshold observed at BES [TJ[28l 

[2^ . The scenario of these states will may be clear in 
future experiments like PANDA at FAIR. 

Beside these, we have also calculated the decay width 
of remaining dimesonic states and listed in Table-IV, 
these states may be identified in the future experimen¬ 
tal measurements. 


5 Conclusion 

In this paper, we are able to calculate the S-wave masses, 
digamma and decay widths of the dimesonic states in 
the light meson sector. Here, we would like to say that 
the dimesonic model becomes more accurate if bind¬ 
ing energy of constituents is small compared to their 
masses. We compared many states which are experi¬ 
mentally observed and predicted theoretically and hav¬ 
ing promising non- qq structure with dimesonic states. 
Calculated results are fairly in good agreement with 
experimental measurements as well as theoretical pre¬ 
dictions. On the basis of mass spectra, decay width and 
in some cases digamma width, we are able to identify 
some dimesonic (meson-antimeson) states. We have cal¬ 
culated the decay widths for experimentally observed 
decay modes for each compared states. 

In our present study, we strongly recommend the states 
/o(980), 5i(1235), /2(1525), /o(1710) as meson-antimeson 
molecules. Furthermore, the states /ii(1380), ao(1450), 
/2(1565), hi (1595) in calculation have been found with 
deviated masses and decay width from experimental 
measurements which rule out their candidature as dime¬ 
sonic molecules. Additionally, some states like /o(1500) 
02 ( 1700 ) and / 2 ( 1810 ) calculated results have little bit 
variation in decay widths, but in some decay modes, 
decay width are fairly near to some experimental re¬ 
sults Some states like 02(1700), /2(1810) 

and /o(1790) have unconfirmed experimental status, 
needing more attentions from experimentalist as well 
as form theoretician. 

Finally, we have predicted the masses, decay widths 
and di-gamma widths of the S-wave dimesonic molecu¬ 
lar states (light sector). Now we would like to employ 
this model to calculate the P-state masses and decay 
properties for light sector as well as for the heavy sec¬ 
tor dimesonic systems. 
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Table 4 Decay width of dimesonic (molecular) states(in MeV). The decay width is calculated for experimentally seen decay 
modes listed in PDG [l]. The other dimesonic states are also tabulated with isospin-spin (I,S). 


States 



Dimesonic decay widths (in 

MeV) 



(in different decay modes) 

total 

PDG yy 

others 

current 





width 



status 








in PDG [J 

/o(980)(,^) 

(49)^^ 

(33)7rp 

( 12 ) 771 . 

94 

50-100 

9itfo Ea 

OK 






108 [2H] 


6i(1235)(^-p) 

(89)^„ 

(43)7r()(, 

- 

132 

142 ±9 

151 ±31 [94] 
113 ± 12 ESI 

OK 


(109)^^ 

(71)^p 

3 

00 

268 

91 ±30 

170 ± 80 1841 

NC! 

ao(1450)(p-) 

(37)7r7j 

(32),,^ 

(25)7r?7' 

94 

265 ± 13 

265 ±30 ESI 
196 ± 10 E7I 

OK 

/o(1500)(pp) 

(50)^^ 

( 29)7777 

(31)kx 

110 

109 ± 7 

108 ± 33 [BS] 

66 ± 10 Ea 

OK 

/'(1525)(^p) 

(25)^^ 

( 23)7777 

( 34)7777 

82 

76± 10 

69 ± 22 [99] 

75 ±4 [Too] 

OK 

/2(1565)(^ij) 

(33)7r7r 

(25)7^k 

( 23 ) 77 ^ 

81 

134 ±8 

113 ±23 Ea 
119 ±24 Ea 

NC! 

hi (1595)^^'7j^ 

(63)u;,7 

- 

- 

63 

384 ± 60 

384 ± 60 UOU 

NC! 


(75)7r77 

(69)7^k 

(37)p7. 

181 

194 ± 40 

151 ± 22 [TSa 
187 ±60 [TSa 

NC! 


(45)x^ 

( 43)7777 

(55)7r7r 

143 

135 ± 7 

1391^2 Esa 

145 ± 8 ESI 

OK 

/2(1810)(^^) 

(75)^^ 

( 62 ) 77,7 

(59) 

196 

197 ± 22 

228^20 Esa 

NC! 

(W)(0,0) 

(756)^^ 

(554)7r7j 

(286)^^ 

1596 

- 

- 

- 

ivv')(o,o) 

(477)^^ 

( 411 ) 777 , 

(355)^^ 

1243 


- 

- 

{v'v')(o,o) 

(841)^^ 

( 709)7777 

(714),,^ 

2324 

- 

- 



(114)^0. 

( 89)7777 

(85)77p 

288 

- 

- 

- 


(74)„^ 

(77).p 

(74) Kx 

225 

- 

- 

- 


(70)„c. 

( 44 ) 7 , 

(70)77P 

184 

- 

- 

- 

{v'p)(i,i) 

(87)^^ 

(87)77p 

(58)xx* 

232 

- 

- 

- 


(175)^^ 

(176)77p 

(139)xx* 

490 

- 

- 

- 

(PP)(0,1) 

(30)7rp 

(30)7rLD 

(15)xx* 

75 

- 

- 

- 

(pp)(r,o) 

(31)^p 

(15)kx* 

( 31 ) 777 . 

77 

- 

- 

- 

(PP)(i,i) 

(29)^p 

( 29 ) 777 . 

(14)xx< 

72 

- 

- 

- 

(PP)(1,2) 

(35),r7r 

(26)p,k 

( 25 ) 777 . 

86 

- 

- 

- 

(tuZJ)(0,0) 

(46),r7r 

(35)p,i^ 

( 33 ) 777 . 

114 

- 

- 

- 

(tuZJ)(0,l) 

(30)7rp 

(15)kx* 

( 29 ) 777 . 

74 

- 

- 

- 


(50)^„ 

(40)p„ 

(36)xx* 

126 

- 

- 

- 

{K*K*)(o,2) 

(60),r7r 

(49)xk 

(47)^7. 

156 

- 

- 

- 

(^‘^*)(l,0) 

(52)^^ 

(43)xx 

(42)^7, 

137 

- 

- 

- 


(48)^.. 

(48)77p 

(35)xx* 

131 

- 

- 

- 

{4>4>)(o,o) 

(21)^^ 

(18)p,k 

(9)x*x" 

48 

- 

- 

- 

{4>4>)(o,i) 

(19)^^ 

(20)77P 

(16)xx< 

55 

- 

- 

- 

{4>4'){o,2) 

(28)T,r 

(25)xx 

(24)^7, 

77 

- 

- 

- 

(P^)(i,r) 

(29)„„ 

(40)7r7r 

(14)xx* 

83 

- 

- 

- 

(P‘^)(l,2) 

(34)7r7r 

(26)p,k 

( 25 ) 777 . 

85 

- 

- 

- 

(<^P)(1,0) 

(51),r7r 

(42)xx 

(41)^7. 

134 

- 

- 

- 

(#)(i,r) 

(47)^„ 

(47).p 

(35)xx* 

129 

- 

- 

- 

(#)(!.2) 

(81)-n-7r 

(67)7ck 

( 64 ) 777 . 

212 

- 

- 

- 

(<()CU)(0,0) 

(53)^^ 

^KK 

( 42 ) 777 . 

139 

- 

- 

- 


(39)7r</, 

(35)ffx* 

(47)^7. 

121 

- 

- 

- 


(NC!- Not Confirmed) 
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